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Abstract 

Spontaneous CP-violation in the minimal supersymmetric stan- 
dard model with a gauge singlet and a cubic superpotential is ex- 
amined. Although the tree-level Higgs potential conserves CP, it is 
shown that with the inclusion of the one-loop top-quark radiative ef- 
fects, CP may be broken spontaneously. The CP-violating minimum 
requires two neutral (/ii,/i2) and one charged (H^) Higgs-bosons to 
be relatively light with m^^ + ~ 100 GeV and mjj± ;$ 110 GeV. 
The electric dipole moment of the electron is in the observable range 
of (I to 3) X 10-2"^ e-cm. 
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There are a number of reasons why it is interesting to study CP violation 
arising from the Higgs sector. Such CP violation could lead to large and 
soon-observable electric dipole moments for the electron and the neutron 
[1]. It may have a role to play in baryogenesis in the early universe [2]. It 
is even possible that this may be the origin of the CP violation observed in 
the Kaon system [3]. And such Higgs sector effects may be readily studied 
at present and future colliders. 

In this Letter we explore CP violation phenomenology in the Higgs sec- 
tor of low-energy supersymmetry, in the simplest model where such effects 
can appear [4]. That model is the minimal supersymmetric standard model 
(MSSM) to which a gauge-singlet Higgs superfield, N, has been added with 
the requirement that the superpotential contain only cubic terms [5]. There 
is no significant CP violation in the parameters of the Higgs effective poten- 
tial (see below), so that any CP violation in the Higgs sector of the model 
must be spontaneous. 

One of the main results we report is that spontaneous CP violation is 
indeed possible in this cubic singlet model. It has been shown by Ramao 
[6] that at tree-level spontaneous CP violation does not occur in this model. 
What allows spontaneous CP violation to occur in our case, in spite of Ra- 
mao's theorem, is the radiative contribution to the effective potential due 
to a heavy top quark [7]. A striking feature of this model is that if CP is 
spontaneously broken then two of the neutral Higgs must be quite light. This 
is actually connected to Ramao's theorem [6]: as the top-quark radiative ef- 
fects are "turned off" the mass-squared of one (or more) of these light Higgs 
goes negative, so that there remains no CP-violating minimum. We show 
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that at best the masses of the two hghtest neutral scalars add up to not much 
more than Mz- In addition, the mass of the charged Higgs is found to be less 
than about 110 GeV. In this scenario the prospects for detecting these Higgs 
particles at LEP and LEP-200 are very good indeed. Moreover, since two of 
the neutral Higgs fields are light, the electron electric dipole moment (EDM) 
comes out in the experimentally accessible range of (10~^^ to 3 x 10"^'') e-cm 
(typically). 

Only the radiative effects of the top-quark will be significant, so that the 
relevant terms in the superpotential are 

W ^XHiH2N+hN^ + htQH2T^ . (1) 

Here //i = ( ) , = ( ) , so that HiH2 = e'^'^H^Hl = (H^^Hl + 
H1H2). The scalar potential for the fields Hi, H2 and N is then given by 
V = Vo + Vtop, where 

vo = l{9i+9imi\'-mr+l9mi\w-\H^H2\') + 

{\HiH2\^ + |iV|'(|i^i|' + 11^21')) + k^\N\^ + Xk{HiH2N*^ + H.c.) + 
XAx{HiH2N + H.c.) + \kAk{N^ + H.c.) + ml\Hi\^ + ml\H2\^ + ml\N\^ 

(2) 

T4op = r^2[ihl\H2? + M%)M{hl\H2\' + M%)/Q'} 

-{hl\H2\Mhl\H2\VQ')}] ■ (3) 

Here = hIHi etc. T^op denotes the leading-log top-quark one-loop 

radiative correction. We have assumed degenerate squarks, M~ = M~ — 
M^g > (174 GeV)^. 



It is easy to see that all the parameters in eq. (2) can be made real by 
field redefinitions, except the ratio r = Ax/Ak- We assume that parameter, 
too, is real. (It would be approximately so in most realistic scenarios [8].) 
To examine the CP-violating minimum, let us define 

D ^ Xk 

E = \A),/\N\ = r\A/\N\ 

F = \kA\N\/\HlHl\ . (4) 

Here and throughout A = A^. Let us further define the phases of the fields 
H^, HI and N as follows: 

avgiHlH^N) = rj 

arg(7V3) ^ ^ . (5) 

There are three terms in eq. (2) that depend on rj and Extremizing the 
potential with respect to these phases gives the CP-violating solution to be 
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Esinr] = — Fsin^ (7) 

These equations can be represented graphically by the triangle in fig. 1. 
Equations (6) and (7) are respectively the law of cosines and the law of sines. 
Denoting the expectation values (VEVs) of the neutral fields by (Hi) — 
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Vi, {H2) = V2 and {N) = x, one can expand the fields around their minimum. 

1 vi 



Hi = vi + 
N ^ x + 



V2\vi\ 

1 V2 

V2\V2\ 
1 X 
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with v"^ = + 1^2 In the absence of CP-violation, $1,2,4 would be 
scalars and $3^5 pseudoscalars. One pseudoscalar gets eaten by the Z. (We 
are working in the unitary gauge.) 

It is straightforward to expand the potential in eq. (2)- (3) about the CP- 
violating extremum to obtain the mass-squared matrix of i = 1, 2.., 5. In 
doing this it is convenient to adopt a common approach of choosing in eq. 
(3) so that at the minimum of Vq the relation V^^^ = is satisfied. (Then 14op 
does not contribute explicitly to the equations for vi, V2 and x, it only corrects 
the (22) entry of the Higgs mass-squared matrix.) The mass-squared matrix 



Xivl + 3rA^V2 

(1 + A)Xivl + SrX'^vl 

—X1V1V2 + (2 - 3r)A^i;i'U2 
Ml, = 

r(4r - 1)A'^ + 36r^X%lvl/x'^ + 4r(6r - l)A^('Ui'y2/a;)cos77 

vfvl/x'^ + 12rXA{viV2/x)cosri 
4:rXA{viV2/x)smr] 
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elements 


are 


Ml 




MI2 




Ml, 




M\2 




Ml, 




Ml, 




Ml 




Ml, 





Ali4 — —rXAv2Cosr] + 2X^xvi — 6rX^V2Vi/x 

= -rXAvxCosr) + 2X^xv2 - 6rX'^vfv2/x 

A^34 = rXAvsinr] 

Ml^ — —3rXAv2smr] 

•^25 — —^rXAvisinr] 

= —3rXAvcosri — 6rX^ViV2v/x (9) 



where Ai = Kfl'i + df) — M'^/v^, A = Aj, and rA = ^4^. Here and henceforth 

Vi,V2 and x stand for the absohitc values of the VEVs. A in the second hne 
of eq. (9) is the top-quark radiative correction defined to be 

The parameter p represents possible non-logarithmic corrections which may 
be significant \i — tji mixing is not negligible [9]. The physical charged 
Higgs boson mass is given by 

m^± = M^, + (3r - 1)AV . (11) 

Since Ai and v — \Jv1-\-V2 — 174 GeV are known, our parameter set is 
{tan/3, X, A, r, A, COST/, A} where tan/3 = V2/V1. (We have used eqs. (4) and 
(6) to eliminate k in favor of the angle cost].) 

Since two of the neutral Higgs bosons {hi and h2) are relatively light in 
this model, their masses and mixings are constrained by the LEP data on 
Higgs searches. In fig. 2 is displayed a plot of the experimentally allowed 
region in the A — cosr] plane for fixed values of the other parameters. The 
first constraint is that hi and /i2 not have been produced in the decay of a 
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real Z [10]. Since there is no significant suppression in the Zh\h2 couphng, 
Z ^ /ii + /i2 should be kinematically forbidden. The bottom (small A) 
boundary of the allowed region in fig. 2 corresponds to m^^+m/ij = Mz- The 
second constraint is that a light Higgs particle not have been produced in the 
decay Z Z*h. \i h = Yl\=i oa^h where X]j = Ij then the cross section 
for the process is approximately proportional to |Q;isin/3 + Q;2COS/3pm^^. The 
condition we have imposed, that gives the upper (large A) boundary of the 
allowed region is [10] mh > (60 Gey)[Q;isin/3 + q;2COS/3|^. In addition, the 
'pseudoscalar' should be heavier than 20 GeV [11], which is always satisfied 
in the fit of fig. 2 (see Table 1). 

The allowed region in fig. 2 corresponds to the fixed parameters having 
the values tan/? — l^r — 0.8,^1 — 3v,a; = 3.8v and A — 3.8. This is 
a relatively large value for A, but is reahzed if the top Yukawa is large, 
near its renormalization-group fixed point. For example, if Mgq — 3 TeV , 
ht = 1.05 (at g2 = (3 TeVy) and p = 1 in eq. (10), we obtain A = 3.8 at 

= mf. (This involves a further running from = to mf assuming 
non-SUSY two-doublet spectrum.) A larger value of ht will be inconsistent 
with perturbative unification (in the desert scenario), so a larger A seems 
to be unhkely, unless the squark masses are heavier than 3 TeV. If A is 
taken to be smaller (with the other parameters fixed at the above values), 
the allowed region rapidly shrinks and vanishes at A ~ 3.1 corresponding to 
Msg ~ 1.5 TeV (with the same ht). Similarly, if tan/? is increased above 1, the 
allowed region shrinks and disappears at tan/3 ~ 1.3. In this scenario of CP- 
violation, then, the parameters A{Msq, ht) and tan/3 are tightly constrained. 
The same is seen from fig. 2 to be true of the parameter A which must lie 
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in a rather narrow range (around 0.2 for r = 0.8). Numerically it is found 
that as r decreases the allowed range of A shifts upward, and vice versa. On 
the other hand, the parameters r,A,x and cosjy can vary over a rather wide 
range. For example, we find solutions for 0.5 < r < 2.7. Much of the above 
behavior can be understood fairly well analytically. 

Consider certain subdeterminants of the 5x5 mass-squared matrix given 
in eq. (9). If the CP-violating extremum is to be a local minimum, this 
matrix must have no negative eigenvalues and its subdeterminants must be 
all non-negative. 



det345 = 3AV^^^x"^(3r - l)sm'^ri{9X%m'^2pv^ + AA^x^^ + 12AAsin2/3cos7/a;i;^} 

(13) 



where detj^ is the determinant of the submatrix corresponding to the ij..k 
rows and columns. Since r > from Ai^-^. cq. (13) implies that (3r — 1) > 0. 
If one sets A = 0, then the right hand side of eq. (12) with (3r — 1) > is 
negative. This is the result of Ramao [6] , that at tree-level spontaneous CP 
violation does not occur in the model. 

Eqs. (12)-(13) allow us to place bounds on A^ if CP is spontaneously 
broken in the Higgs sector. 



One also sees from eq. (12) and (13) that the hght Higgs masses are made 
largest by making sin2/3 large, that is, tan/? ~ 1 is preferred. As a conse- 
quence of eq. (14), the charged Higgs boson mass is predicted to be less than 



det 



1235 





(14) 
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about 110 GeV (see cq. 11). 

More insight can be obtained by considering the small A limit, which is 

justified by eq. (14) and by our numerical results. As A — > 0, the mixing 

between the $1,2,3 block and $4,5 block goes away. If one block diagonalizes 

to obtain the 3x3 block to order A^ one finds 

Xivj -i;it;2(Ai + 2(3r - 1)A2) 0\ 

A<2($i,$2,$3) = I -^;i^;2(Ai + 2(3r - 1)A2) X^vUl + A) 

0/ 



+ CAi 



^ v|cOS^77 l'ii;2COS^77 — OT2Sin7^COS7y \ 



t;it;2COS^77 v\cos'^r] —vvismrjcosrj 

\ —vv2sm'r]C0sr] —vvismrjcosr] v'^sin^r) 



j (15) 



where C = [4r(3r - l)/(4r - 1)] (A^/Ai). From the eigenvalues of A4^(<l>i, $2, $3) 
a non-trivial lower bound on A^ can be derived as follows. Label the eigen- 
values to be Mf < Mjj < Mjjj. It is possible by considering the eigenvalues 
as a function of C (defined above) to show that M| < CM| = CXiv"^. (This 
can be done by plotting the eigenvalues of M.'^{C)/Mz versus C and consid- 
ering how often they cross the line y = C which can be determined from the 
equation dct{7Vl^(C)/M| — CI} = 0. This is easily seen to be quadratic (not 
cubic) in C and to have two real roots, one of which is C = 0.) From eq. (15) 
one sees that for small A^ and A > cot^/3— 1 the largest entry of A4^ is greater 
than (1 + A)Xiv^ = (1 + A)M|sin2/3. Therefore Mfjj > (1 + A)M|sin2/3. 
The trace of implies that Mj + Mjj + M]jj = M| ( 1 + Asin^/J + C) . Thus 
M| + Mjj <Ml{C + cos^ 13). This together with Mj < CM^ imply that if 
C < cos2/5, then {Mj + Mu)/Mz < + cos/3, while if C > cos^/?, then 
{Mi + Mil) /Mz < V2C + 2cos^ 13. In order to have Mi + Mii>Mz it then 
must be that C > (1 — cos/3) ^. Using the definition of C one gets a lower 
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bound on A^, which combining with eq. (14) gives 

^^^^^Ai(l - cosPf < (3r - 1)A' < ^XiiVl + A - 1) . (16) 

For the values of r = 0.8, A = 3.8, tan/3 = 1 that are used in fig. 2 one has 
0.13 < A < 0.41 which is seen to be satisfied well by the numerical results. 
One can also derive from eq. (16) a better lower limit on r than r > 1/3. 

A remark about the small A hmit is in order. Prom eq. (4) and fig.l it 
can be seen that as A — > 0, x/A r/k. {k in fig. 2 is about 0.63, the fixed 
point value.) Terms have been dropped in deriving eq. (15) that are down by 
0{Xx/A). For large r, these terms become important, which tend to decrease 
the fight eigenvalues oi Ai^. This is why no realistic solution exists for large 
r (r > 2.7). 

In fig. 2 we have also plotted contours of the electric dipole moment of the 
electron. These have been computed using the results of Leigh, Paban and Xu 
(LPX) [1] . Only the graphs with the same topology as the top quark graphs 
(figs. 1 and 2 of LPX), and only the contributions from the lightest two 
neutral Higgs have been included. This should be good to about 20%. The 
typical values of de are in the range (10~^^ to 3 x 10~^^) e-cm, interesting for 
future atomic experiments. In Table 1 are presented the masses and mixings 
of the lightest Higgs for various A and cost] in fig. 2 from which we have 
computed de- 

It is interesting to ask if the Higgs sector CP-violation alone can account 
for e and e' in the K meson system. Pomarol [3] has shown that this may 
indeed be possible in a more generalized version of the SUSY singlet model. 
Since the effects of spontaneous CP violation felt in the fermion sector in our 
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model is identical, we expect the results of Ref. 3 to hold in our case as well. 

As for the neutron edm (ci„), although it is difficult to predict its value 
precisely due to strong interaction uncertainties, we expect dn to be near 
the experimental limit, d„ ~ 10~^^ e-cm. Indeed, if the only source of CP- 
violation is in the Higgs sector, then the usual fine-tuning of the phase of the 
gluino mass term can be avoided. This phase will be zero at tree-level, and 
will arise only at one-loop, giving experimentally consistent value for dn- 

So far we have said nothing about whether the CP-violating minimum 
that we have been studying is also the absolute minimum. For fixed |i?2| 
and |A^| it is easy to sec that V at the CP-violating extremal angles (eq. 6) 
is lower than at the CP conserving values r], C = 0, tt. Of course the true 
CP-violating and CP-conserving minima will not be at the same values of 
Hi, H2 and N. However, it can be shown analytically that for small A the CP- 
violating minimum is the lowest [9] . We have verified this result numerically 
for some reasonable values of the parameters. 
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Figure Captions 



Fig. 1. This triangle shows the relation between the extremal values of 
the phase angles r] and ( and the magnitudes of the fields and the parameters 
of the Higgs potential, li D,E and F are not in such a relation that this 
triangle can be drawn, then there is no CP-violating extremum. 

Fig. 2. A plot of the allowed region (inside the closed curves) of the 
(A — COST]) plane, for tan/? = l,r = 0.8, A = 3v,x = 3.8v and A = 3.8. 
The lower boundary of the region corresponds to m/u + m/jj — Mz- The 
upper boundary is the constraint from non-observation oi Z ^ Z* + h. The 
contours inside the region represent the electron EDM in units of 10~^^ e-cm. 

Table 1. Some selected points in the (A — cosrj) plane for the same param- 
eter values as in fig. 2. m^^ and ruh^ are the masses (in GeV) of the hghtest 

two neutral Higgs hi and /i2- The Rij are defined by hi = Z)f=i Rij^j- The 
electron EDM de is in units of 10~^^e-cm. 
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